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ABSTRACT 

We present the first collective evidence that Fermi-detected jets of high kinetic power (ikin) are 
dominated by inverse Compton emission from upscattered external photons. Using a sample with 
a broad range in orientation angle, including radio galaxies and blazars, we find that very high 
power sources (Lkin > 10^^'^ erg s~^) show a significant increase in the ratio of inverse Compton 
to synchrotron power (Compton dominance) with decreasing orientation angle, as measured by the 
radio core dominance and confirmed by the distribution of superluminal speeds. This increase is 
consistent with beaming expectations for external Compton (EC) emission, but not for synchrotron- 
self Compton (SSC) emission. For the lowest power jets (Lkin < lO^'^'^ erg s~^), no trend between 
Compton and radio core dominance is found, consistent with SSC. Importantly, the EC trend is 
not seen for moderately high power flat spectrum radio quasars with strong external photon fields. 
Coupled with the evidence that jet power is linked to the jet speed (Kharb ct al. 2010), this finding 
suggests that external photon fields become the dominant source of seed photons in the jet comoving 
frame only for the faster and therefore more powerful jets. 
Subject headings: galaxies: active — quasars: general — radiation mechanisms: non-thermal 



1. INTRODUCTION 

The relativistic jets of radio-loud active galactic nu- 
clei (AGN) are copious gamma-ray emitters, as was first 
discovered by EGRET ([Hartman et al.|[l992D . and con- 
firmed by the Fermi large area telescope (LAT), which 
has associated over 800 s ources with radio-loud A GN in 
the second catalog f2LAC: lAckermann et al.l2011[) . Most 
of these are blazars, seen with the jet axis along the line 
of sight, though Fermi has also detected the jets of sev- 
eral r adio galaxies (RG; Abdo et al. 2 0101 : iKataoka et al.l 
120 lit) , which are mis aligned blaz ars under the standard 
unification scheme (Urrv fc Padovani 1995iV 

The lower-energy peak in the jet spectrum is well- 
understood as synchrotron emission from relativistic 
electrons in the jet. The high-energy component, peak- 
ing from X-ray to TeV energies, is attributed to pho- 
tons upscattered by the same relativistic electrons to 
higher en ergies via the inverse Compton (IC) pro- 
cess (see iBottcheil 120071 for a review). These pho- 
tons could arise from the jet synchrotro n emission (syn- 
chrotron sclf-Compton emission, SSC; IMaraschi et al.l 
[1992; Marschcr & Travis 1996) o r from externa l sources 
such as the accretion disk (Dormer et al. 1992), broad- 
line r egion (BLR: .Sikora et al...l994.'). or rnolccu lar torus 
(;MT: IBlazeiowski et al.ll200a : ISikora et al.ll20Q9D . i.e. ex- 
ternal Compton (EC) emission. Identifying the IC emis- 
sion mechanism is a diagnostic for the location of the 
gamma-ray emittin g region, a currently open issue (e.g. 
lAgudo et al.ll20li[ ). However, the spectral energy dis- 
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tributions (SEDs) of individual sources are rarely suf- 
ficiently constraining of th e IC mechanism d ue to the 
number of free parameters (jSikora et al.lll997l) . 

The beaming pattern (how apparent luminosity 
changes with orientation) is diff'erent for EC and SSC 
emission. Thus, a collective study using a sample of 
sources at different orientations can be used to identify 
the gamma-ray emission mechanism. In this letter we 
discuss the effect of jet power and orientation on the ob- 
served IC power in view of the recent suggestion of a 
dichotomy in radio- loud AGN (iJH). We then show that 
the collective beaming pattern for a sub-set of high-power 
sources appears to support EC models rather than SSC 
for this population (^SJ, and discuss the implications for 
the gamma-ray emission region. 

2. A DICHOTOMY IN RADIO-LOUD AGN 

2.1. The Synchrotron Plane 

Using a large sample of jets, we recently found ev- 
idence that radio-loud AGN form two populations in 
the plane of synchrotron peak luminosity (Lp^syn) ver- 
sus peak frequency (f^p.syn; Meyer ct al. 2011, hereafter 
Mil). A population of 'weak' jets consists of sources 
with low -/jp,syn which appear to be most aligned at high 
f^p.syn (^10^^ Hz) and trace out a shallow track on the 
plane as they become less aligned (dropping more in 
i^p.synthan Lp,syn)- A Separate population of 'strong' 
jets with higher ip_sy„ and t'p.syn ^ 10^^ Hz appear to 
drop rapidly in luminosity with misalignment (as mea- 
sured by radio core dominance). Importantly, the weak 
jets exhibit jet kinetic powers below lO'*''''^ erg s~^ while 
all sources with Lkin > lO'*^'^ erg s~^ are on the strong- 
jet branch. The weak/strong divide in the synchrotron 
plane, which we associate with the morphological di- 
chotomy in Fanaroff-Riley (FR) type 1 and 2 RG, may 
also be mapped to a critical transition in accretion effi- 
ciency (.Ghiscllini fc Cclotti.,2001; ,Gliiscllini et al..,2009al: 
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Fig. 1. — This plot of total LAT-band luminosity versus gamma- 
ray energy index a (approximating Lp.ic versus i^p,ic) is the high- 
energy analog of the synchrotron plane discussed in Mil. Blazars 
have been divided based on the luminosity and location of their 
synchrotron peaks into strong-jet and weak-jet blazars as described 
in the text. Those without !/p,aynand Lp^syn are 'unclassified'. 
FR2 RG appear to have similar !-'p,syn (similar a) to the strong- 
jet blazars while FRI RG have much lower values of fp.syn (higher 
a) relative to the weak-jet blazars with which they are associated. 
The FRI would be detectable by Fermi even at much harder a (see 
lAckermann et al.ii2011l . Figure 15). The region at upper right re- 
mains empty despite re cent claims of high-power, high-peak sources 
(shown as black stars. iPadovaniet al..,2012, ). matching what has 
been found in the synchrotron plane. 

iGeorganopoulos et al.l[20Tll ). 

2.2. The Blazar Envelope at High Energies 

In light of the dichotomy discussed above, we present in 
Figure[l]the total LAT-band luminositjO (a proxy for the 
IC peak luminosity, ip,ic)i versus the LAT-band energy 
index a (a proxy for peak frequency, Vp.ic) for the entire 
Fermi-detected sample of radio-loud AGN of known red- 
shift (data taken from the 2LAC). Sources with higher 
z/pjc will have harder (smaller) values of a, thus in the 
figure v-pjc increases to the right. 

An empty region at upper right appears analogous to 
that seen in the synchrotron plane. IGiommi et al.l (|2012[ ) 
have suggested that this is a selection effect due to a lack 
of redshifts for high-frcqucncy-peaked, high-luminosity 
sources, and lPadovani et al.i (2012) discuss four such can- 
didates. However, these candidates do not cross into the 
upper-right region in Figure [1] (black stars). These inter- 
esting sources appear to exhibit t'p^syn at or above ~ 10^^ 
Hz; however, the SED sampling is sparse, and it is dif- 
ficult to rule out an alternative explanation for the soft 
X-ray spectra such as an extra emission componen t (see 
e.g., the case of BL Lacertae: iRaiteri et al.l[2010( ). Im- 
portantly, for the several hundred Fermi sources lacking 
identifications or redshifts, nearly all (>99%) have soft 

^ iiOOMcV-lOOGcV = '^T^d^S^i / {1 + z)^~°' for power-law sources, 
where 5-, is the catalog integrated energy flux and a = F— 1 using 
the published photon index F. For sources fit with Log-Parabolic 
spectra, the integrated luminosity is calculated from the 2LAC 
fitted values with a k-correction. 



spectra (a > 1) or would require redshifts > 2 or higher 
to cross into this region, making it unlikely to be empty 
due to selection effects. 

The blazars in Figure [1] have been divided in an ap- 
proximate way based on their synchrotron spectra into 
strong and weak jets. We take the former to be those 
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well as flat spectrum radio quasar (FSRQ) type sources 
outside this area. For the weak jets, we take all sources 
outside this area, minus any FSRQ. While very rare cases 
of FRl-like (e.g., weak) FSRQs exist, their occurrence is 
negligible in the brig ht fermt-detected p opulation con- 
sidered here (see e.g.. iKharb et al.l 120101 and references 
therein). The eleven Fermz- detected RG are also shown 
(data from lAbdo et al.ll2010D . The FRI in Figure [1] have 
much lower J^pjc (larger a) compared with the weak-jet 
sources, while there is little difference between a for FR2 
and strong jets, tentatively matching the different mis- 
alignment paths for these populations that were found in 
the synchrotron plane. 

2.3. The Importance of Jet Kinetic Power 

In Mil, we showed that Lkinis an important param- 
eter in classifying RL AGN. We have selected from the 
2LAC a subset of 152 blazars (as well as the detected 
radio galaxies) with estimates of the 300 MHz isotropic 
lobe emission, which is scaled to estimate Lkin as in Mil. 
For 46 sources with adequate coverage of the high-energy 
SED, we estimated Lpjcfrom a two-sided parabolic fit 
to the 2LAC data in combinatio n with X-ray data take n 
from NEEQ and/or Swift /^KT (|Cusumano et alllMoh . 
For an additional 106 sources, ip ic is scaled from the 
rest-frame luminosity at 1 GeV estimated from the 2LAC 
fitted spectrum: log Fp,ic = log -^iGcV + ??(«— 1)^, where 
77=1.5 if a > 1, 77 = 5 otherwise. The difference between 
the two estimates is typically less than 0.2 decades in log 
L. 

We use the radio core dominance i?cc= log (^corc/^oxt) 
as a tracer of the orientation angle, where the core lu- 
minosity is measured at 1.4 GHz and the extended at 
300 MHz (see Mil). As discussed in Mil, the absolute 
normalization between orientation angle and i?co depends 
on Lkinj with i?ce decreasing with increasing orientation 
angle for a given Lkin- 

In Figure [2] we show the 145 sources which have a 
known redshift. ipjc is positively correlated with both 
i?ccand ikin- The OLS bisector fits to the combined 
blazar-RG sample shown have slopes of 1.1, 1.5, 1.4, 1.3, 
1.5 (±0.1) from lowest to highest bin in Lkin- The corre- 
lation between ip jc and i?ce in each group is significant 
and positive (R values from 0.37 to 0.89); the apparent 
gamma-ray output of a blazar is therefore a strong func- 
tion of both Lkin and the orientation angle. 

2.4. Emission Mechanisms and Their Beaming Patterns 

For both sy nchrot ron and SSG, the beaming pattern is 
L = L'JP"*"" (jDerm er 1995), where L assumes isotropic 
emission in the galaxy frame, L' is the solid-angle in- 
tegrated luminosity in the comoving jet frame, 5 is the 
Doppler factor, and a the spectral index. For the expo- 
nent p+a, p=3 for a 'moving b lob', or p=2 for a st ation- 
ary feature in a continuous jet (jLind fc Blandfordl[l985) . 

® NASA Extragalactic Database: |http://ned.ipac.caltech.edu/| 
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TABLE 1 
Predicted Correlation Slope Values 



index values 



ssc 



EC 



Psyn P'y ^7 ^cc ^p JXci 



1.1-1.3 
1.6-2 
1.6-2 
1.2-1.5 



1.7-2 
2.4-3 
2.4-3 
1.7-2.5 



0.6-0.7 
0.8-1 
1.2-1.5 
0.8-1 



with the possibihty of different values for different emit- 
ting regions. For th e EC cas e, the pattern is Lbc = 
L'^c6P+^+^" (|Dermed[l995HGeorganoDOulos et al.l[200l . 

The slope of Lpjc versus Rcc depends on the gamma- 
ray emission process. From the synchrotron beaming, 

Rce = l0g(Lcore/iext) = (Pr + Ctr) log(5 -|- Ci , (1) 

where the factors Cn depend on the unbeamed luminosity 
and jet power, but do not affect the slope. We take a = 
1 for the peaks, and eliminating S we obtain: 



log Lpoak 



Re 



C2, 



(2) 



where h = I + p^ (SSC) or b = 3 + p^ (EC). For ref- 
erence we list in Table [1] the expected slopes for various 
cases of p values assuming typical values of a^ (0—0.5). 
The slopes in Figure [2] appear to be more consistent with 
SSC under a p^, Pr=3 or p^, Pr=2 scenario. However, 
the lack of i^ermi-detected misaligned FR2 sources (par- 
ticularly, with _Rco < —0.5) is likely to affect the slope 
for the two highest bins in Lkinin Figure [2] we show for 
reference the upper limits to ipjcQ for four RG in the 
Mil sample. We take the rightmost of these points as 
the most constraining: including it in the fit for the log 
^kin = 45—45.5 bin increases the slope to 1.5±0.2, and 
for the highest bin to 1.7±0.3. These slopes, and our as- 
sumption that the Lorentz factor (F) of the gamma-ray 
emitting plasma is similar to that in radio, are discussed 
below. 

3. THE COMPTON DOMINANCE IN FERMI BLAZARS 

3.1. EC in powerful FSRQ 

For sources with a high-energy component dominated 
by SSC, the Compton dominance (i?p ) should remain 
constant over all orientations (no correlation between 
i?p and i?cc expected) . For the EC case, 

Rp = log{LEc/Ls) = {p-, - Psyn + 2) logS + C3. (3) 

Using Equation [1] we have the relation 

P7 Psyn I ^ 



Rp — 



Pr 



Re 



C4. 



(4) 



As shown in Table [TJ we expect a correlation with a slope 
from 0.6—1.5 under the simplest assumptions for sources 
emitting gamma-rays by EC. 

Figure [3] shows Rp versus Rcc for the sample discussed 
in Section 12.31 Sources are divided into broad bins of 



'' From the flux limit (> 100 MeV): log /loo = -8.3 + 
(5/6)(«— 1.2), with average a =1.3 for our sample and Lpjc scaled 
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Fig. 2. — Estimated IC peak luminosity (X/p,ic) versus radio core 
dominance (i?co) for 145 blazars and ten RG (circled) detected 
by Fermi. The correlation between ipjc ^nd Rcc becomes clear 
when binning by Lkin ■ OLS bisector fits to sources grouped by 
L]jin have slopes 1.1, 1.5, 1.4, 1.3, 1.5 (±1; Pearson's Rvalues 0.81, 
0.89, 0.82, 0.37, 0.44) beginning with the lowest bin. For an EC 
emission model, these slopes are predicted to be higher (> 1.7); 
however, for the high-power sources, the lack of Fermi detections 
for misaligned (i.e., low Rcc) sources will affect the apparent slope. 

Lkin- These bins include a mix of sources in terms of op- 
tical type (FSRQ or BL Lac); however all the very high 
power (VHP, Lkin > 10^^'^) are 'strong-jet' types in our 
classification due to their Lkin and position in the syn- 
chrotron plane, including many apparent BL Lacs, which 
likely suffer dilution of their broad lines by th e jet emis- 
sion (iGeorganopoulos fc Marsched l T9 98: Ghis ellini et alJ 
12011'). For the moderate-power sources, the division into 
strong/weak roughly follows the FSRQ/BL Lac divide. 

No trend between Rp and i?ce is evident for blazars 
overall. However, when VHP sources are selected, a pos- 
itive correlation emerges. We also show the upper lim- 
its on i?p for VHP sources from the sample of Mil not 
detected by Ferm^. VHP sources from the 2LAC and 
Mil with upper limits on LcxtU sue also shown; because 
of the binning on Lkin which is scaled from Loxtj we also 
show as connected gray points the maximum Rce such 
that Lkin > 10''^-^ erg s~^ The OLS bisector fit through 
all VHP points (including upper/lower limits) gives a 
slope of l.lzbO.l; this is a lower limit since most of the 
Rp upper limits are on the lower half of the correlation. A 
sub-sample of sources with Rce values contemporaneous 
to Fermi, calculated from the average 15 GH43 flux over 
the time-frame of the 2LAC as m easured by the Owen s 
Valley Radio Observatory ( OVRO. [Richards et al.ll201lD . 
gives an identical slope. 

The non-zero slope suggests EC emission in VHP 
sources; however, when taken together the slopes in Fig- 

' See previous footnote, with index a estimated from the relation 



for detected sources: a=—0.2xur, 
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^ Estimated at 300 MHz from the lowest-frequency SED point 
with spectral index a = 1.2. 



from LiGeV f Section 12.3 



^° scaled to 1.4 GHz by -1.2 in log L 
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Fig. 3. — The Compton dominance (-Rp) versus radio core dom- 
inance (-Rce). A positive correlation is seen for the most powerful 
sources (shaded), as predicted for EC emission. The slope of the 
OLS bisector fit to the shaded points is l.litO.l, though the pres- 
ence of the upper limits on Rp and lower limits on Rcc (range shown 
by black lines) makes this a lower limit on the real slope. Weak 
jets (black points) as well as moderately powerful FSRQ (open 
triangles) show no trend, as expected for simple SSC. 



ures [3] and [3] are inconsistent with a simple EC scenario, 
as seen from Table[TJ We have assumed that S is the same 
for all energies; if the radio-emitting plasma is slower, 
then the expected slopes above will increase from the 
quoted ranges for both figures. Assuming truly higher 
slopes for the VHP sources in Figure [51 concordance can 
be achieved with a general p=3 model with a slower F in 
the radio. 

3.2. The Test of Superluminal Motions 

As shown by iLister fc Marscheil (ll_997) , a substantial 
number of sources are expected to be seen at orientation 
angles smaller than that oi = 1/F which maximizes 
the superluminal speed Papp- We therefore expect the 
highest /3app values at Rp {Rcc ) values less than the max- 
imum. Indeed if the highest -Rp (-Rcc ) sources are most 
aligned (within 1/F), their /3app values should be rela- 
tively small (noting that these may be under-represented 
as their /3app are difficult to measure). In Figure |4] (upper 
panel) the VHP sources are plotted as in Figure [31 col- 
ored according to /3app . The highest /3app are observed 
in the middle of the correlation, as expected. This is 
also seen in the lower panel, where we have estimated 
the maximum /3app in bins of width 0.25 in -Rcc, running 
parallel to the correlation line (3 

3.3. Why is EC Only Apparent in Powerful Jets? 

What is apparently depicted in Figure [3] is that 
strong jets at high -L^in sue EC dominated, while lower 



The estimator for /3a 



is derived from the order statis- 



tic Y'„=max(/3app) for a bin of size n using the estimator for 
the unknown upper bound of a uniform distribution /9app,max = 
Y„{n + l)/n with variance /3^pp,max/("-^ +2n) ('e.g. lCasella|[2051) . 
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Fig. 4. — Upper: The same plot as Figure [S] for high-power 
sources (ikin > lO'*^'^ erg s~^) with measurements of maximum 
parent iet s peed (/3app). Val u es are taken from [Jp rstad ct al] 



t spe! 

J20Q11, [20051); IK ellermann et al.l | [200l) : ILister et a O (2009), and 
;Kharb et al. (201^). The appearance of low /3app sources with high 
core and Compton dominance supports the interpretation that 
these sources are aligned very close to the line of sight. Lower: 



Sources are binned along the shown correlation (running perpen- 
dicular to the dashed line in units of 0.25 of -Rcoas illustrated by 
the two dotted lines), and the maximum /3app is estimated. The 
increase in /3app, max from lower .Rcc /iJp towards a maximum, fol- 
lowed by a decrease towards the highest values is qualitatively as 
expected if -Rce and -Rp increase with alignment. 

-Lkin strong jets are not. However, the typical error in log 
-Lkinof 0.7 is a significant factor, as it may be the case 
that the restriction Lkin>10*^'^ erg s~^ is simply that 
which is high enough to avoid any contamination with 
SSC sources. Thus the true boundary Lkin for the SSC to 
EC transition is probably lower than 10""^-^ (indeed, a few 
FSRQ of moderate -Lkin are co-spatial in Figure [31 with 
the VHP sources). However, our findings imply that (a) 
VHP sources are dominated by EC, and (b) many lower 
Lkin strong jets appear at a similar (low) -Rp regardless of 
-Rcc, suggesting that EC is not important for some part 
of the population. 

Which emission mechanism dominates is related to 
how fast the flow is: assuming the gamma-ray emission of 
strong jets is inside the BLR or MT, the comoving energy 
density of the external photon field is F^f/cxt- The syn- 
chrotron energy density is Us = Ls/{4:7rc^t'^^^T^), where 
ivar is the observed variability timescale in hours, Lg the 
synchrotron luminosity, and we have assumed S—T. EC 
will dominate over SSC provided F^C/cxt > Us, or for F 
greater than a transition value 



Ftr = 16.2 



is/10 



l47 



(ivar/6)'(C/ext/10-4)^ 



(5) 



Typical estimates of C/oxt in the BLR and MT differ by 
~ 100 (t/ext,BLR ~ 2 6 X 10"^ erg cm~^ an d /7ext,MT « 
2.6 X IQ-^ erg cm-3: lGhisellini et al.ll2009b( ). In the MT 
we have Ftr,MT ~ 14.4, higher by a factor of 100^/^ = 1.8 
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than in the BLR, where rtr,BLR '^ 8.1. 

However, VLBI studies indic ate that F > 10 for most 
FSRQ (e.g. lJorstadet aLll2005l) : if the GeV emission site 
is in the BLR, it is difficult to explain the lack of EC sig- 
nature for some strong jets. It has also been found that 
sources with higher -/ ^kin produce on the average faster su- 
perluminal motions (jKharb et al.ll2010() . The connection 
between L^in and F, and an emitting region in the MT, 
can then explain a transition to EC at high Lkin : Strong 
jets at lower Lkin are also slower and as long as F < Ftr 
they are SSC emitters, exhibiting a Compton dominance 
independent of radio core dominance. As Lkin increases, 
F also increases and once F > Ftr, the photon field of 
the MT begins to dominate, producing a Compton dom- 
inance that increases with increasing radio core domi- 
nance, as seen in Figure [31 If we adopt longer ivar, both 
Ftr, BLR and Ftr,MT dccreasc, requiring that essentially 
all strong jets are EC emitters, which our data does not 
support . On the other han d, with a shorter ivar ~1 hr 
(see e.g. lFoschini et aLll2011|) we obtain Ftr — 12.7 for the 
BLR and Ftr = 22.5 for the MT, making either location 
plausible. 

4. CONCLUSIONS 

The gamma-ray luminosities of Lermi-detected radio- 
loud AGN have been shown to depend strongly on both 
Lkin and the orientation angle. We find the first collective 
evidence for external Compton emission in high-power 
jets (Lkin > 10^^-^ erg s~^): as can be seen in Figure [31 



the Compton dominance of these sources increases with 
radio core dominance, a measure of orientation. This re- 
quires that the beaming pattern of the gamma-ray emis- 
sion is more focused than that of synchrotron, as is the 
case for EC scattering. A confirmation that Rp increases 
with decreasing orientation angle comes from the fact 
that the apparent superluminal speeds are observed to 
increase along the correlation to a maximum at moder- 
ate values of Rp, before decreasing towards the highest 
values (Figure[3]), as expected since we anticipate a signif- 
icant number of sources at angles smaller than the maxi- 
mum superluminal speed angle 1/F. The fact that strong 
jets of lower Lkin do not show an increase of Compton 
dominance with alignment suggests that their gamma- 
ray emission is due t o SSC. Because m ore powerful jets 
appear to be faster (jKharb et al.l 120101 ) . a transition to 
EC at high Lkin can be explained by a transition Lorentz 
factor above which the external photons dominate. 
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